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Abstract 

In adaptating to perpetual darkness, cave species gradually lose eyes and body pigmentation and evolve alternatives for 
exploring their environments. Although troglodyte features evolved independently many times in cavefish, we do not yet 
know whether independent evolution of these characters involves common genetic mechanisms. Surface-dwelling and 
many cave-dwelling species make the freshwater teleost genus Sinocyclocheilus an excellent model for studying the 
evolution of adaptations to life in constant darkness. We compared the mature retinal histology of surface and cave 
species in Sinocyclocheilus and found that adult cavefish showed a reduction in the number and length of photoreceptor 
cells. To identify genes and genetic pathways that evolved in constant darkness, we used RNA-seq to compare eyes of 
surface and cave species. De novo transcriptome assemblies were developed for both species, and contigs were annotated 
with gene ontology. Results from cave-dwelling Sinocyclocheilus revealed reduced transcription of phototransduction and 
other genes important for retinal function. In contrast to the blind Mexican tetra cavefish Astyanax mexicanus, our results 
on morphologies and gene expression suggest that evolved retinal reduction in cave-dwelling Sinocyclocheilus occurs 
in a lens-independent fashion by the reduced proliferation and down regulation of transcriptional factors shown to have 
direct roles in retinal development and maintenance, including cone-rod homeobox (crx) and Wnt pathway members. 
These results show that the independent evolution of retinal degeneration in cavefish can occur by different develop- 
mental genetic mechanisms. 
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Introduction 

Cave animals evolve distinct troglomorphic characters in per- 
petual darkness, most commonly, the reduction of eyesight, 
pigmentation, pineal organ, and scales, countered by aug- 
mented chemoreceptors, mechanoreceptors, and lipid stor- 
age (Poulson and White 1969; Jeffery 2001, 2009). Although 
these traits evolve in common across many cave taxa, an 
open question remains: Is the evolution of troglomorphic 
traits due to the common and repeated evolution of ortho- 
logous genetic and developmental mechanisms? The most 
extensive investigation of the evolution of troglodyte traits 
in cavefish comes from investigations of the Mexican tetra 
Astyanax mexicanus. In cave Astyanax, lens apoptosis associ- 
ated with overexpression of sonic hedgehog (shh) induces 
eye degeneration (Yamamoto et al. 2004; Alunni et al. 2007; 
Borowsky 2008). In Mexican tetra cavefish, mutations in oca2 
and mclr appear to be responsible for the loss of eye and 
body pigmentation and the "brown" body color morph, re- 
spectively, (Protas et al. 2006; Gross et al. 2009). Although 
Astyanax is an excellent model for understanding the molec- 
ular mechanisms of regressive evolution, research on a variety 
of cave-animal models is necessary to understand whether 
independent evolutionary lineages utilize related molecular 



genetic mechanisms, or if evolution can achieve similar results 
by different mechanisms. 

The investigation of retinal degeneration in cavefish is not 
only important to improve our understanding of the evolu- 
tionary mechanisms of development but also because it rep- 
resents a significant evolutionary mutant model for human 
disease (Albertson et al. 2009). Age-related macular degener- 
ation (AMD) is the primary cause of blindness in old age and 
is predicted to affect approximately 3 million people in the 
United States by 2020 (Klein et al. 1997; Gehrs et al. 2006). 
AMD involves the loss of the retinal pigmented epithelium 
and photoreceptors, a phenotype remarkably similar to that 
found in some cavefishes (McCauley et al. 2004; Gross et al. 
2009). Although advances have been made in treating the 
vascular aspects of AMD, little progress has been made for 
other forms of the disease (Ambati et al. 2003). The investi- 
gation of the genetic mechanisms underlying evolutionary 
models of retinal degeneration may provide clues that may 
lead to potential new therapies for AMD. 

China has a rich cavefish fauna with currently 101 valid 
species that belong to a single order (Cypriniformes) compris- 
ing only three families; Cyprinidae (including zebrafish, gold- 
fish, carp, and 56% of all Chinese cavefishes), Cobitidae (true 



© The Author 2013. Published by Oxford University Press. 
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http:// 
creativecommons.org/licenses/by-nc/3.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, 
provided the original work is properly cited. For commercial re-use, please contact journals.permissions@oup.com 



Open Access 



Mol. Biol. Evol 30(7):1 527-1 543 doi:10.1093/molbev/mst079 Advance Access publication April 23, 2013 



1527 



Mengetal. • doi:10.1093/molbev/msc079 



MBE 



loaches, 2% of cavefishes), and Balitoridae (river loaches, 42% 
of cavefishes) (Zhao et al. 2011). Southwestern China is one 
of the largest cave-rich karst geomorphologic regions in 
the world (Huang et al. 2008). Different karst types provide 
suitable conditions to support a high diversity of cavefish 
populations. The freshwater teleost genus Sinocyclocheilus 
(Cypriniformes: Cyprinidae) is endemic to the karst region 
of the east Yungui Plateau and northwest Guangxi in south- 
western China. The genus contains over 55 known species, 
including many surface-dwelling species and at least 10 cave- 
dwelling species with different degrees of eye degeneration 
and melanin loss (Xiao et al. 2005; Chen, Zhang, et al. 2009). 
Situations such as this with gradients of visual function pos- 
sess advantages for analyzing the mechanisms of retinal de- 
generation (Tobler et al. 2010). Mitochondrial DNA 
sequences show that Sinocyclocheilus species lie within 
the Cyprinion-Onychostoma lineage, related to rock carp 
(Procypris rabaudi) and common carp (Cyprinus carpio) 
(Wu et al. 2010) (supplementary fig. S1, Supplementary 
Material online). Common carp, goldfish, and related species 
have an extra round of whole genome duplication compared 
with zebrafish (Danio rerio) (Wang et al. 2012). The genus 
Sinocyclocheilus likely shares that tetraploid origin and fish of 
this genus have 96 chromosomes, about twice that of most 
teleosts (Xiao et al. 2002). High species diversity and pheno- 
typic variation make Sinocyclocheilus an ideal model for 
the investigation of evolutionary questions. We report here 
integrated transcriptome-wide investigations of changes in 
gene expression patterns related to the evolution of eye 
degeneration in Sinocyclocheilus cavefish. 

To evaluate morphological evolution of cavefish eyes, we 
used antibody markers to detect specific retinal cell types in 
tissue taken from adult blind cavefish (Sinocyclocheilus 
anophthalmus, the "blind goldenline barbel") and compared 
these with a closely related surface species (S. angustiporus, 
the "small gill opening goldenline barbel"). To understand 
evolved changes in gene expression, we utilized next- 
generation sequencing technologies (Wang et al. 2009). 
The Sinocyclocheilus transcriptome, which we assembled de 
novo from lllumina sequencing reads, provides a valuable 
resource for future studies, including the facilitation of studies 
on other species in this genus and the construction of a 
gene reference for annotating an eventual Sinocyclocheilus 
genome assembly. In this study, we identified eye genes 
that were differentially expressed between surface and cave 
species, and found enriched signaling pathways and disease- 
associated genes among differentially expressed genes. Results 
showed that evolved retinal degeneration in cave-dwelling 
Sinocyclocheilus is not associated with lens disappearance, 
in contrast to A. mexicanus, but is more likely due to the 
down regulation of transcription factors shown to have 
direct roles in retinal development and maintenance, includ- 
ing cone-rod homeobox (crx) and Wnt pathway members. 
These results suggest that convergent troglodyte characters 
can evolve by different genetic mechanisms in different 
fish taxa. Candidate genes identified by transcriptome-wide 
analysis may give clues to the mechanisms controlling 
the evolution of eye degeneration in cavefish and provide 



an alternative evolutionary medical model (Albertson et al. 
2009) for retinal degeneration that often accompanies 
old age. 

Results 

Sinocyclocheilus anophthalmus Cavefish Evolved 
Decreased Eye Size 

Sinocyclocheilus anophthalmus lives in constant darkness and, 
as its name suggests, lacks external eyes; it has a duckbill-like 
snout, sinking frontal bone, and long barbels and pectoral fins 
compared with surface fish. Sinocyclocheilus anophthalmus 
do not show a detectable response to natural light, even 
after stimulation with a laser pointer (Zhao and Zhang 
2009) (fig. 1A). Consistent with the small internal eyes of 
adult S. anophthalmus, we discovered that its optic lobe is 
reduced in size compared with surface species (Meng FW, 
Braasch I, Zhang CG, Postlethwait JH, unpublished 
data). Sinocyclocheilus angustiporus is a close relative of 
S. anophthalmus, which has normal eyes, and although 
occasionally entering caves, usually occupies surface waters 
(fig. 1B) (Zhao and Zhang 2009). Sinocyclocheilus anophthal- 
mus retains small eyes that are buried deeply within adipose 
tissue and are covered with skin (fig. 1A and A'). The exam- 
ination of different developmental stages was not possible 
in this study due to the unavailability of cavefish embryos. 
But, to assess the extent of eye reduction in adult fish, we 
dissected the fish and measured eye diameters. Results 
showed that the internal eye of S. anophthalmus cavefish 
was about a third the size of the external eyes of adult surface 
fish (fig. 1C and D), although the size of retinal cell bodies, 
as assayed by both nuclear and cell-specific labeling, is com- 
parable between the two species. Hematoxylin and eosin 
(H&E)-stained sections showed that the eyes of cavefish 
retained the basic vertebrate eye structure, with fully formed 
lens, cornea, iris and neural retina (fig. 2A and B). Although 
the cavefish neural retina retained its layered cellular 
organization, it was notably thinner than that of surface 
fish, with fewer cells in all three nuclear layers (fig. 2C 
and D). The morphological changes of photoreceptors in 
Sinocyclocheilus cavefish appeared to be more substantial 
than that of inner retinal neurons in H&E-stained sections. 

The vertebrate retina contains six major types of neurons 
organized as laminar structures, including rod and cone pho- 
toreceptors, as well as horizontal, bipolar, amacrine, and gan- 
glion cells; these cell types form parallel microcircuits that 
integrate and process visual signals. To assay the distribution 
and numbers of different specific retinal cell types in eyes of 
cavefish and surface fish, we used the monoclonal antibodies 
Zpr1 which labels red/green double cones and Zpr3, a marker 
for rod outer segments and a subset of cone outer segments 
(Yin et al. 2012). To visualize cell types in the inner retina, we 
used an antibody against protein kinase C, which labels ON 
bipolar cells in the inner nuclear layer (INL), and an antibody 
against glutamine synthetase, which labels Muller glia. We 
found that cones, rods and bipolar cells were all present in 
cavefish eyes, with overall cell body sizes indistinguishable 
from those of surface fish. However, their numbers were 
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Fig. 1. Cavefish phenotypes. Cavefish Sinocyclocheilus anophthalmus (A) and surface fish S. angustiporus (B). Cavefish have no external eyes (A), 
but they have internal eyes (A r and C, left) that are much smaller than the external eyes of their surface fish relatives (C, right). (D) Comparison of 
eye diameter normalized with total fish body length including the caudal fin. Values expressed as mean ± SD, ** P < 0.01. Scale bar in (A and B): 1 cm; 
(C): 2 mm. 



significantly reduced (fig. 2E-J). In cavefish eyes, cones (fig. 2E 
and F) and rods (fig. 2C and H) were 49% and 59% shorter 
than in surface fish eyes (P < 0.001, cave n = 6, surface n = 7, 
t test), respectively. The presence of photoreceptor outer seg- 
ments was confirmed by H&E staining (fig. 2C and D) and by 
zpr-3 labeling (fig. 2C and H), although the lateral organiza- 
tion of the outer segments was disrupted in the cavefish 
retinas. In the inner retina, Muller glial cells were reduced in 
cavefish compared with surface fish (fig. 2K-L). Counts of 
nuclei indicated that the cell density of the outer nuclear 
layer (ONL), INL, and ganglion cell layer (CCL) were all signif- 
icantly decreased in cavefish relative to surface fish (fig. 2A/I). 
Cell number reduction of the ONL, which contains the 
photoreceptor nuclei, was more substantial than that of the 
INL and GCL (fig. 2A/I). We conclude that the evolution of 
Sinocyclocheilus cavefish involved a reduction in the number 
of photoreceptors, bipolar cells, and Muller glia, as well 
as substantial morphological abnormalities in photoreceptor 
structure. 

Cell Proliferation and Apoptosis in the Cavefish Retina 
We next tested whether cell proliferation in cavefish eyes was 
affected during adaptive evolution, using an antibody 
directed against phosphorylated histone-3 (PHH3), an 
M-phase marker. In the eyes of surface fish (fig. 3A and C), 
mitotic cells were detected throughout the retina, with the 
highest levels of proliferation localized to the ciliary marginal 
zone (CMZ), in which retinal progenitors are located 
(Raymond et al. 2006) (1.91 ± 0.33 positive cells per linear 
mm of retina, n - 5). In cavefish retina, however, the overall 
number of proliferating cells was reduced (0.79 ± 0.31 positive 



cells per linear mm of retina, n = E>; P< 0.001) with signifi- 
cantly fewer mitotiocally active cells noted in the CMZ as 
compared with the surface species (fig. 3B and D). We con- 
clude that the cavefish eye is not replenished as the fish grows 
as happens in surface fish. 

To analyze apoptosis in adult surface and cave species, we 
performed antibody labeling with anti-active Caspase-3, an 
early marker of apoptosis (Dai and Krantz 1999; Ryu et al. 
2005). No apoptosis was detected in the retina of either spe- 
cies (fig. 3E and F). As a positive control, we investigated 
induced cell death in zebrafish retina; results showed that 
Caspase-3 antibody detected cell apoptosis induced by ultra- 
violet (UV) light in 3dpf zebrafish (fig. 3G). Although apopto- 
sis was not detected in adult cavefish, we cannot exclude the 
possibility that excess apoptosis occurred in the embryonic or 
larval stages of cavefish eyes relative to surface fish eyes. 

Sequencing, De Novo Assembly, and Annotation of 
the Sinocyclocheilus Reference Transcriptome 
To help understand differences between surface and cavefish 
eyes on a genetic level, we profiled gene expression. We iso- 
lated total RNAs from eyes and brains of S. anophthalmus and 
S. anophthalmus (two individuals for each species), and con- 
structed four cDNA libraries: cavefish eyes, cavefish brain, 
surface fish eyes, and surface fish brain. Two brain cDNA 
libraries of both species were sequenced using the lllumina 
Genome Analyzer llx platform to produce 80 nucleotide (nt) 
paired-end reads (SRA run accession numbers: surface fish, 
SRR788094; cavefish, SRR788095), and two eye cDNA libraries 
of both species using lllumina HiSeq 2000 with 100-nt paired 
end reads (surface fish, SRR788096; cavefish, SRR788097). 
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Fig. 2. Eye degeneration in cavefish. H&E stained sections of adult surface fish eye (A, C) and cavefish eye (B, D). Representative double-labeling in the 
retina of surface (E, G, /) and cavefish (F, H,j) eyes; cell nuclei are dye-labeled red with ToPro3; retinal cell types are labeled green (anti-mouse Alexa Fluor- 
488). Cones (Zpr1; F) and rods (Zpr3; H) of cavefish were shorter than those of surface fish (E, G). Dotted brackets show the length of cones and rods. 
ON bipolar cells were labeled with antibodies against protein kinase C (PKC, green; I, J). MCiller glial cells were labeled with antibodies against glutamine 
synthetase (GS, green; K, L). Cavefish had significantly fewer nuclei in the ONL, INL, and GCL (A/I) normalized for eye size. Values expressed as 
mean ± standard deviation (SD); *P < 0.05; **P < 0.01. C, cornea; I, iris; NR, neural retina; RPE, retinal pigmented epithelium; IS, inner segment; OS, outer 
segment; BC, biplolar cell; MG, MCiller glial cells. Scale bar in (A, B): 1 mm; (C-L): 50 urn. 



After trimming adapters and removing low quality reads, a 
total of 132,682,597 reads, containing 9,179,015,231 nt, were 
retained (supplementary table S1, Supplementary Material 
online). De novo assembly on surface and cave species 
using Trinity (Grabherr et al. 2011) assembled 66,459,580 
reads into contigs with 300 nt minimum length from two 
species leaving 66,223,017 unassembled reads (singletons). 
The contig file for cavefish contained 152,464 contigs and 
that for surface fish contained 156,118 contigs (supplemen- 
tary table S1, Supplementary Material online). To generate a 
common and nonredundant Sinocyclocheilus transcriptome 
database for further RNA-seq analyses, we merged the 
two contig files using Cap3 (Huang and Madan 1999). 
The combined Sinocyclocheilus transcriptome contained 
56,074 contigs, which were on average longer than from the 

1530 



single-species assemblies (supplementary table S1 and fig. S2, 
Supplementary Material online). 

Next, we compared the combined Sinocyclocheilus tran- 
scriptome with the UniGene records for zebrafish (D. rerio), 
medaka (Oryzia latipes), and three-spined stickleback 
(Gasterosteus aculeatus). Blastn similarity searches (cut-off 
E value of 10 -3 ) showed that 30,590 Sinocyclocheilus contigs 
matched 19,724 zebrafish UniGenes, covering 37.5% of 
the zebrafish UniGene records. Medaka and stickleback 
databases yielded fewer matches than did the zebrafish tran- 
scriptome: 1,794 Sinocyclocheilus contigs matched 1,531 
medaka UniGenes (6.5%), and 1,853 Sinocyclocheilus contigs 
matched 1,554 stickleback UniGenes (8.3%). Two factors likely 
contribute to finding more UniGene hits in zebrafish than 
in stickleback and medaka: first, a phylogenetic tree based on 
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Fig. 3. Cavefish have reduced retinal cell proliferation. (A-D) Retinas from adult surface fish (A, C) and cavefish (B, D) were stained with anti-PHH3, a 
marker of the M-phase of the cell cycle. The CMZ lies at the anterior edge of the peripheral retina. Proliferating cells were located mainly in the CMZ in 
surface fish retinas (A). PHH3 positive cells were not detected in the CMZ of cavefish (B). PHH3 positive cell numbers per linear millimeter of retina were 
markedly reduced in non-CMZ of cavefish retinas (C: 0.79 ± 0.31, n = 5) relative to that of surface fish (D: 1.91 ± 0.33, n = 5, P < 0.001). (E-G) Caspase-3 
staining was performed to analyze apoptosis in adult surface (£) and cave (F) eyes. Apoptosis was not apparent in the retina of either surface or cave 
species, although in positive controls, Caspase-3 antibody did detect apoptosis induced by UV in 3dpf zebrafish eyes (G). OS, outer segment. Scale bar in 
(A-F): 50 Jim. 



de novo and downloaded sequences of z/'c family member-1 
(zicl) and C protein-coupled receptor-85 (gpr8S) shows that 
Sinocyclocheilus is more closely related to D. rerio than any 
other species in the tree (fig. 4), and second, the zebrafish 
UniGene database (92,836 entries) is substantially larger than 
that of medaka (22,128 entries) and stickleback (16,737 
entries). 

For annotation, contigs were searched against the Nr data- 
base using blastx. We found that 38,017 contigs (67.8%) had 



at least one hit (E-value cut-off of 10 -6 ), while 18,057 contigs 
(32.2%) did not align to any known protein in the Nr data- 
base; many of these likely represent 3 r and 5' untranslated 
regions (UTRs). The 38,017 contigs that hit known proteins 
grouped into 20,062 UniGenes. Most contigs (31,407, 82.6%) 
aligned best to D. rerio UniGenes. To globally characterize the 
transcriptome, we annotated Sinocyclocheilus contigs based 
on sequence homology to functionally annotated sequences 
in other species. We identified 30,721 contigs (54.8% of 
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Fig. 4. Phylogenetic position of Sinocyclocheilus. A phylogenetic tree generated by Maximum Likelihood analysis of the combined nucleotide sequences 
of the CGD paraiogs of zid and gpr85, which have a single copy in human and zebrafish. Analysis confirms that the Sinocyclocheilus lineage experienced 
a genome duplication after its divergence from zebrafish and shows that zebrafish is the closest species to Sinocyclocheilus with a sequenced genome. 
Numbers at the nodes denote the bootstrap percentages of 100 replicates. Sequences for all taxa can be found in supplementary table S5, 
Supplementary Material online. 



Sinocyclocheilus contigs) that had affinity to at least one Gene 
Ontology (GO) term. In addition, of the 38,017 contigs with 
good hits to the Nr database, 7,296 were classified as un- 
known, predicted, hypothetical, or uncharacterized proteins 
(fig. 5). 

We assigned contigs to broad functional categories and 
found that 30,721 contigs were associated with 124 generic 
GO-Slim annotations; specifically, these included 35 catego- 
ries at the Cellular Component level, 39 categories at the 
Molecular Function level and 50 categories at the Biological 
Process level. 

Differential Gene Expression in Surface Fish and 
Cavefish Eyes 

To investigate the evolution of the cavefish eye transcrip- 
tome, we compared surface fish and cavefish transcriptomes. 
Eye samples of both species were sequenced with lllumina 
HiSeq 2000, and resulting reads were trimmed for adapter and 
quality, aligned to the Sinocyclocheilus transcriptome, 




No Blast Hits No Mapping No Annotation Annotated TotaJ 



Fig. 5. Annotation of the Sinocyclocheilus transcriptome. Blast2GO was 
used to annotate our de novo assembled transcriptome using param- 
eters described in Materials and Methods. A total of 30,721 contigs 
(54.8% of total contigs) were annotated; 13% of contigs had no GO 
term but at least one blast hit in the Nr database. Many contigs that 
failed to annotate may represent 5 r - or ^-UTRs. 
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counted, and normalized to the total mapped reads. For the 
cavefish sample, 22,793,525 reads (49.31%) mapped to 
Sinocyclocheilus transcriptome contigs and an additional 
412,525 unmapped reads aligned to zebrafish UniGenes. 
(Unmapped reads were those that had not assembled into 
Sinocyclocheilus contigs.) A total of 23,206,050 reads from 
eye samples, corresponding to 49.9% of all high quality 
reads in cavefish, mapped to 52,915 contigs (94.4%) in the 
Sinocyclocheilus transcriptome. For the surface fish sample, 
29,584,957 reads (48.98%) mapped to the Sinocyclocheilus 
transcriptome, whereas an additional 562,988 reads mapped 
to zebrafish UniGenes; a total of 30,147,945 surface fish reads 
(50.3%) aligned to 52,656 Sinocyclocheilus reference contigs 
(93.9%). To analyze gene expression profiles, we converted 
read numbers mapped onto UniGenes to RPKM (reads per 
kilobase of exon per million mapped sequence reads) 
(Mortazavi et al. 2008). 

Among 9,649 unique genes whose RPKM values were over 
5 in the Sinocyclocheilus eye transcriptome, we identified 
1,658 genes as differentially regulated, defined as greater 
than a 2-fold change between cavefish and its related surface 
species with a P-value less than 0.05 (log 2 -fold changes 
[log 2 FC] > 1, P value < 0.05). Of differentially expressed 
genes, 1,189 (71.7%) had decreased expression and 469 
(28.3%) had increased expression in cavefish eyes (fig. 6A). 
Supplementary table S2, Supplementary Material online, lists 
genes according to their change in expression level (log 2 FC). 
To validate the differential expression of genes identified by 
RNA-seq, we used real-time polymerase chain reaction (PCR) 
with primers derived from the de novo assembled 
Sinocyclocheilus transcriptome to quantify the mRNA levels 



of several genes that were selected from different change 
levels and different expression levels. Real-time PCR data 
validated the RNA-seq results (fig. 6B; supplementary fig. S3, 
Supplementary Material online), showing that differential 
gene expression in cavefish is reliable and that the de novo 
assembled Sinocyclocheilus transcriptome can be used for 
RNA-seq based gene expression analysis. 

Analysis showed that the RNA-seq results confirmed 
conclusions from histological investigations. We found that 
photoreceptor genes were downregulated (rhodopsin, 
guanine nucleotide binding protein [G protein], alpha trans- 
ducing activity polypeptide 7 [gnatl], gnat2, and opsin 7 [cone 
pigments], long-wavesensitive, 7 and short-wave-sensitive 
2 [opnllwl and opnlsw2] with log 2 FC = -1.5803, -1.7425, 
-1.7428, -2.1126, and -1.4817, respectively), as was the 
Muller glia gene glutamine synthetase (log 2 FC= —1.3598). 

Identification of Enriched Pathways and Diseases 
among Differential Expressed Genes 
Several GO categories were over- or underrepresented in the 
differentially expressed genes in each of the three main GO 
categories (supplementary fig. S4, Supplementary Material 
online). For example, in the biological processes category, 
signal transduction (GO: 0007165) and protein modification 
process (GO: 0006464) were overrepresented in the downreg- 
ulated gene group, likely reflecting the low number of 
photoreceptors present in cavefish eyes (fig. 2). The protein- 
aceous extracellular matrix (GO: 0005578) and extracellular 
space categories (GO: 0005615) were overrepresented in the 
upregulated gene group in the cellular component category, 
likely reflecting the larger spaces between cells observed in 




^^^^^^^^^^^^■^^^■^■^^^^■H Fold change of Cave to Surface fish 

Log2 (RPKM of surface fish eye) 

Fig. 6. (A) Comparison of two eye RNA-seq determinations measured in reads per kilobase of exon per million mapped sequence reads (RPKM). Lateral 
green lines bracket genes with 2-fold differences in expression; these genes were considered to be differentially expressed in this study. A total of 469 
genes were upregulated and 1,189 were downregulated in cavefish relative to surface fish. The purple line is the regression line (R 2 = 0.706). White circles 
indicate genes verified by real-time PCR (B) using RNA isolated from eyes of each species. The expression levels of atpSal, crxa, crygsl, hsp47, nr2e3, otx2, 
pdeSaa, rhodopsin, and sox2 were quantified and normalized to beta-actinl . Relative expression values are mean ± SD of at least three independent 
experiments. 
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our histological studies. Seven of 13 genes encoded by the 
mitochondrial genome were downregulated in cavefish eyes 
(mt~col, 8.1 x; mt-nd6, 6.8 x; mt-ndS, 5.5 x; mt<yb, 5.2 x; mt- 
nd4l, 2.7 x; mt-co2, 2.5 x; and mt-nd4, 2.2 x), suggesting de- 
creases in rates of energy metabolism in cavefish eyes. 

To further explore altered regulation of signal pathways 
during cavefish evolution, we used KOBAS 2.0 (Xie et al. 201 1 ) 
to identify enriched pathways and disease-associated genes 
among the differentially expressed genes (supplementary 
table S3, Supplementary Material online). Twenty-two path- 
ways and 13 diseases had significant P values less than 0.05 
after FDR correction in the downregulated gene group and 18 
pathways and 23 diseases were significantly enriched in the 
upregulated group. Although identified pathways had some 
overlap among different databases within KOBAS 2.0 (Xie 
et al. 2011), we found that the most significantly enriched 
pathways were relevant to visual phototransduction in the 
downregulated group, including phototransduction (KEGG), 
visual signal transduction: cones and rods (PID Curated), syn- 
aptic transmission (Reactome), and heterotrimeric G-protein 
signaling pathway-rod outer segment phototransduction 
(PANTHER). Most of the enriched diseases were eye pathol- 
ogies related to the retina. For example, retinitis pigmentosa 
(Fun DO), retinal disease (Fun DO), cone-rod dystrophy 
(GAD), achromatopsia (GAD), and bradyopsia (OMIM). In 
the upregulated group, enriched pathways and diseases 
were diverse, involving cytokine-cytokine receptor interac- 
tion (KEGG), muscle contraction (Reactome), cell adhesion 
(KEGG), pertussis (KEGG), cardiac diseases (KEGG), Wnt sig- 
naling pathway (PANTHER), and immune system diseases 
(KEGG). These transcriptomic results are as would be pre- 
dicted from the histological analyses showing a strong reduc- 
tion in photoreceptors. 

Expression of Retinal Development and Maintenance 
Genes in Cavefish 

To identify potential regulatory differences in retinal cell 
development and maintenance genes in adult 
Sinocyclocheilus cavefish, we queried genes already shown 
to be involved in fish retinal development. During photore- 
ceptor development, six key transcription factors affect rod 
and cone differentiation and maintenance, including RAR- 
related orphan receptor beta (rorb), orthodenticle homolog-2 
(otx2), crx (for which zebrafish has two co-orthologs of 
the tetrapod gene, currently called crx and otxS but which 
should be called crxa and crxb[otxS], respectively [Suda et al. 
2009]), neural retina leucine zipper (nrl), nuclear receptor 
subfamily 2 group E member 3 (nr2e3, also called retina- 
specific nuclear receptor) and thyroid hormone receptor 
beta (thrb) (Swaroop et al. 2010). 

Vertebrates possess five groups of light-sensitive opsin 
proteins that mediate the conversion of photons into an 
electrochemical signal (Fain et al. 2010). Opsin expression is 
regulated by photoreceptor-enriched factors, including crx, 
retinal homeobox (rx), nr2e3, nrl, and thrb, or ubiquitously 
expressed factors, such as ataxin 1-like 3 (atxnll3), K(lysine)' 
acetyltransferase 2B (kat2b), E1A binding protein p300 a 



(ep300a), and ep300b (Peng and Chen 2007; Hennig et al. 
2008). Our RNA-seq experiments showed that cone opsins 
and rhodopsin were significantly downregulated in cavefish 
eyes and real-time PCR confirmed these expression differ- 
ences (fig. 6; supplementary fig. S3, Supplementary Material 
online). Immunostaining studies of rhodopsin in tissue 
sections provided results expected from real-time PCR inves- 
tigations: rhodopsin was distributed throughout the photo- 
receptor layer in surface fish eyes, but localized in a thinner 
layer in cavefish eyes (fig. 2E and F). In some cavefish, changes 
in opsin gene expression represent an early evolutionary step 
toward retinal degeneration (Tobler et al. 2010). Several tran- 
scription factors involved in opsin expression were downreg- 
ulated more than 2-fold: rx2 (-1.09 log 2 FC), nr2e3 (-1.24), 
crxa (-1.45), and crxb (otxS, -1.17). In contrast, rxl (-0.03 
log 2 FC), nrl (0.32), thrb (—0.32), and ubiquitously expressed 
transcriptional factors, atxnll3 (0.03), kat2b (-0.08), ep300a 
(—0.15), and ep300b (0.35), were expressed at similar levels in 
both species (fig. 7; supplementary table S2, Supplementary 
Material online). 

Experiments knocking down crxa gene expression in 
zebrafish had shown that many genes in the visual transduc- 
tion cascade were downregulated, including rhodopsin, cone 
opsins, interphotoreceptor retinoid-binding protein (irbp), 
otxS(crxb), gnatl, gnat2, and guanylate cyclase activator 
lA(gucala) (Shen and Raymond 2004). Similarly, in 
S. anophthalmus cavefish, the expression levels of a series of 
phototransduction genes, including genes reduced in crxa 
morpholino treated zebrafish, were also decreased (fig. 7; 
supplementary table S2, Supplementary Material online). 
These data support the hypothesis that the reduction of 
crxa expression that we observed in cavefish eyes gives rise 
to the downregulation of phototransduction genes in cave- 
fish eyes. We also designed primers to detect the expression 
level of the two different crxa paralogs in Sinocyclocheilus that 
derived from the carp genome duplication (CGD), crxaa 
and crxab, with real-time PCR (fig. 8A). Results showed that 
both crxaa and crxab were downregulated in cavefish relative 
to surface fish, but the ratio between crxaa and crxab was 
similar in each species (fig. 8B). 

Comparison between Sinocyclocheilus and Astyanax 
Cavefish 

Eye degeneration in A mexicanus cavefish involves lens cell 
apoptosis induced by shh overexpression (Yamamoto et al. 
2004). In contrast, the lens did not show obvious morpholog- 
ical changes in adult S. anophthalmus cavefish, (fig. 2A and B). 
Lens measurements in histological sections of Sinocyclocheilus 
eyes revealed no substantial difference between the species 
regarding the relative areas of the lens normalized by the area 
of eye sections (35.8 ± 8.4% in surface fish; 35.2 ± 9.0% in cave- 
fish). Although the expression of shh also did not change 
significantly in the eyes of adult S. anophthalmus relative to 
the eyes of surface species (log 2 FC= —0.01), it is the expres- 
sion of shh from surrounding tissues during embryonic 
development, not shh expression in the eye itself, that is 
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Fig. 7. Proposed model for the evolution of eye degeneration in the cavefish Sinocyclocheilus anophthalmias. (A) Surface fish. (B) Cavefish. In normal 
vertebrate eyes (A), otx2 upregulates crx, which controls the differentiation of photoreceptor precursors; a number of other factors, including nr2e3, nrl, 
thrb, and crx itself, help to regulate the differentiation of precursor cells into mature rods and cones that express differentiation products, including 
gnatl, gnatl, gucala, and irbp (Hennig et al. 2008; Swaroop et al. 2010). In cavefish eyes (B), transcriptomic data showed downregulation of crx, which 
would lead to decreased photoreceptor cell differentiation, fewer and malformed photoreceptors, and diminished expression of a series of photore- 
ceptor-specific genes (rod opsins and cone opsins, otx5, gnatl , gnatl, gucal a, and irbp). Because otxl expression was unchanged in cavefish eyes, another 
regulatory factor may lie upstream of crx, and/or evolved changes in the promotors of crx paralogs in cavefish may dampen the response of crx to 
normal otxl signal. Because fish can regenerate photoreceptors from the CMZ and from MCiller glia cells, this pathway likely works throughout life and is 
thus involved in retinal maintenance. Downregulated genes in cavefish are indicated by green letters. Upregulated genes included several immune- 
related genes, which may represent response to inflammation associated with eye degeneration. 
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Fig. 8. (A) Our de novo assembled Sinocyclocheilus transcriptome identified two co-orthologs, crxaa and crxab, of the single zebrafish crxa gene; these 
arose in the tetraploid origin of Sinocyclocheilus. The crxaa and crxab genes have several identifying nucleotide variations. To distinguish the expression 
level of different crxa paralogs, we designed primers to detect crxaa and crxab with real-time PCR using RNA isolated from eyes of each species. (B) The 
expression levels of crxa paralogs were quantified and normalized to beta-actin-1 . Relative expression values are mean ± SD of at least three independent 
experiments. Results showed that both crxaa and crxab were downregulated in cavefish, but that there was no obvious difference between crxaa and 
crxab within each species. 



important for eye development in Astyanax (Yamamoto et al. 
2004; Pottin et al. 2011). 

Next, we compared our RNA-seq results from adult 
Sinocyclocheilus cavefish and surface fish with expression 



studies that used a zebrafish microarray to query gene expres- 
sion from whole bodies of 3-day-old larvae of A mexicanus 
cavefish and surface fish (Strickler and Jeffery 2009). 
Comparisons revealed that 28 of the 67 differentially 
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expressed genes in A. mexicanus larval bodies were repre- 
sented in our RNA-seq adult eye transcriptome gene set 
(fig. 9). Some of the 39 missing genes may have been 
expressed in the many cell types other than the eye found 
in the whole body of Astyanax larvae. Of the 28 genes present 
in both data sets, only rhodopsin, cytokine inducible SH2- 
containing protein (cish), gnatl, and gnatl, had greater than 
2-fold change between cave and surface species in our 
experiments (boxed in fig. 9). Several crystallin genes were 
expressed in the epithelial cells of the lens, and downregulated 
in Astyanax cavefish larvae (Strickler and Jeffery 2009). The 
anti-apoptotic activity of aA-crystallin was implicated in 
Astyanax eye degeneration (Strickler et al. 2007; Strickler 
and Jeffery 2009), but the ctA<rystallin gene was not present 
in the Sinocyclocheilus transcriptome, presumably because it 
may be expressed only in early development (Strickler et al. 
2007). The expression levels of several crystallins, including 
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srpr 
ptpro 



Fig. 9. Cluster analysis of genes expressed differentially between cavefish 
and surface fish comparing Astyanax and Sinocyclocheilus. The Astyanax 
transcriptome data came from the cross-species hybridization of 
Astyanax larval whole body cDNA to a microarray of zebrafish 
sequences (Strickler and Jeffery 2009), which identified 67 genes with 
differential expression. Comparisons between Astyanax and 
Sinocydocheilus revealed that, of 28 genes in the Astyanax data set 
that were also detected in the Sinocyclocheilus RNA-seq data, only 
four were differentially expressed in both species of cavefish. These 
four genes (rhodopsin, cytokine inducible SH2~containing protein [cish], 
gnatl, and gnatl), are indicated by blue circles. Color intensity repre- 
sents upregulation and down regulation, respectively. 



y/\A2~crystallins and /3B~crystallins, were actually higher in 
Sinocyclocheilus cavefish than in surface fish (supplementary 
table S2, Supplementary Material online), providing further 
evidence that lens degeneration is unlikely to be driving 
retinal degeneration in Sinocyclocheilus cavefish. 

Discussion 

Retinal Defects in Cave Sinocyclocheilus 
We analyzed morphological and transcriptomic differences 
in the eyes of closely related cavefish and surface fish in the 
genus Sinocyclocheilus to suggest genes, genetic pathways, 
and developmental mechanisms that evolved during the 
adaptation of S. anophthalmus cavefish to constant darkness. 
Results showed that cavefish eyes were imbedded in fat under 
the skin, and were reduced in overall size. Retinal layers were 
normally patterned in cavefish, but were reduced in thickness 
and cell density, especially in the photoreceptor layer. 
Although retinal cell bodies were about the same size in 
cavefish and surface fish, the apical-basal length of photore- 
ceptor cells was significantly shorter in cavefish compared 
with surface fish (fig. 2E-H), and photoreceptor outer seg- 
ments were disorganized. In addition, retinal cell proliferation 
was greatly reduced in cavefish eyes relative to surface fish 
eyes, and gene expression patterns were substantially altered 
in cave and surface species, especially the expression of genes 
previously shown to regulate retinal development. 

Phylogenetic analyses that combined morphology and mi- 
tochondrial gene sequences divided genus Sinocyclocheilus 
into four clades (Zhao and Zhang 2009). The cave and surface 
species studied here, S. anophthalmus and S. angustiporus, 
both belong to Clade tingi, which diverged from other 
clades before the Second Act of the Tibetan Orogeny about 
2.6 million years ago. Some species within this clade 
inhabit caves and others occupy surface streams belonging 
to the same water system, located only a few kilometers 
from each other. S. anophthalmus is one of the few 
troglophile species in Clade tingi, and likely separated from 
other species in this clade about 1 million years ago (Zhao and 
Zhang 2009). 

Histological investigations showed that, although cavefish 
eyes retained the gross retinal architecture of vertebrate 
eyes with lenses of approximately normal size relative 
to the size of the eyeball, they had thin retinas, short photo- 
receptor cells, and more sparsely populated retinal cell 
layers compared with surface fish. Although we were not 
able to obtain cavefish at early developmental stages 
for this study, we detected several classes of differentiated 
retinal cells organized into distinct retinal cell layers in adult 
cavefish, indicating that embryonic eye induction, differenti- 
ation, and patterning signals are intact during cavefish 
development. 

In contrast to mammals, the eyes of adult individuals 
of cyprinids and other teleost fish grow continuously and 
can regenerate in response to retinal damage (Otteson and 
Hitchcock 2003). In mature zebrafish, new photoreceptors 
arise from two sources, the CMZ at the periphery of the 
retina and the Miiller glia cells in the INL (Raymond et al. 
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2006). Although no obvious apoptosis was found in cavefish 
(fig. 3£ and F), cell proliferation rates were reduced in cavefish 
photoreceptors (fig. 3B and D) compared with surface fish, 
most notably in the CMZ region (fig. 3A and C). This obser- 
vation suggests the hypothesis that the reduction in cavefish 
eye size results from changes in the genetic regulation of 
proliferative signals in the retina during in postembryonic 
stages. The morphology of the adult cavefish eye, in which 
individual retinal cell body size is comparable with that of 
surface fish, but cell height, cell number, and overall eye 
size are dramatically reduced, provides further evidence 
for evolutionary changes specifically affecting growth and 
proliferation factors. 

Data from our transcriptomics experiments further 
support conclusions inferred from morphological analyses. 
Among 9,649 unique genes, we identified 1,189 that were 
down regulated and 469 that were upregulated in cavefish 
eyes relative to surface fish eyes. The major pathways and 
disease-associated genes that were downregulated in eyes 
of cavefish versus surface fish involved visual transduction 
of the retina (including protein modification, visual signal 
transduction, G-protein signaling, respiratory electron trans- 
port, and synaptic transmission). These results were consis- 
tent with our observations on eye morphology showing 
abnormal retinal cell morphology. Genetic factors associated 
with retinal disease in humans, including retinitis pigmentosa, 
cone-rod dystrophy, achromatopsia, and macular degenera- 
tion, were also notably downregulated. 

We conclude that developing cavefish do not lack embry- 
onic eye induction signals, because all normal histological cell 
types are present in adult cavefish eyes and these layers enjoy 
nearly normal spatial relationships to each other. Although 
cavefish eyes possess the same cell types as surface fish, cell 
morphologies are altered and cell numbers are diminished. 
We conclude that the proliferation and/or maintenance 
of retinal cells, especially photoreceptor cells, is defective in 
cavefish. This conclusion suggests that the degeneration of 
cavefish eyes results from changes in genetic regulation 
during late stages of retinal development and/or in the main- 
tenance and postembryonic growth of eyes rather than in 
events specifying early stages of eye development. 

Gene Pathways Altered in the Cavefish Eye 
Transcriptome 

Several gene pathways were altered in cavefish eyes relative to 
those of surface fish. Results showed that wnt2 (4.03 x), lefl 
(3.16x) and, importantly, several Wnt antagonists, including 
axin2 (2.06x), naked cuticle homolog 7 (nkdl, 2.1 3 x), sfrpla 
(1.55x), and sfrp2 (2.64x), were upregulated in cavefish eyes. 
The Wnt signaling pathway directs growth and differentiation 
in retina and other tissues (Van Raay and Vetter 2004; Clevers 
2006). Analysis of model systems suggests that the Wnt 
pathway plays an anti-apoptotic role in the retina (Sundin 
et al. 2005; Yi et al. 2007; Chen, Hu, et al. 2009). The Wnt 
antagonists Sfrpl, 2, 3, and 5 are upregulated in human retinas 
affected by retinitis pigmentosa, an inherited form of photo- 
receptor degeneration that generally begins in the peripheral 



retina with rod cell dysfunction. Upregulation of Sfrp2 reduces 
Wnt signaling, leading to increased cell death in the retina 
(Jones et al. 2000a, 2000b; Yi et al. 2007). Although we did not 
detect apoptosis in the retinas of adult cavefish, we cannot 
rule out the possibility that cell death at an earlier time 
point might play a role in reducing retinal size in adult 
Sinocyclocheilus cavefish. The observed upregulation of sfrp2 
presents the possibility that reduced Wnt signaling in cavefish 
retinas may lead to photoreceptor degeneration at an earlier 
stage of life. The significance of the upregulation of wnt2 in 
cavefish eyes in the face of the upregulation of Wnt inhibitors 
warrants further study. 

The photoreceptor-specifk nuclear receptor nr2e3 was 
downregulated in Sinocyclocheilus cavefish eyes. Mutations 
of the human NR2E3 gene are responsible for several retinop- 
athies, including clumped pigmentary retinal degeneration, 
which involves night blindness, rudimentary rod function, 
and hyperfunction of the "blue" S-cones (Schorderet and 
Escher 2009). NR2E3 activates transcription of several rod- 
specific genes, including rhodopsin (Cheng et al. 2004), and 
suppresses cone genes (Chen et al. 2005; Hennig et al. 2008; 
Swaroop et al. 2010). The down regulation of nr2e3 in 
Sinocyclocheilus cavefish could contribute to the reduction 
in opsin expression we observed in RNA-seq analysis. 

Results showed increased expression of z/c7, zic2a, and 
zic2b in cavefish eyes by 5.7x, 2.2x, and 7.5x, respectively, 
whereas zic3 expression remained unchanged. The embryonic 
mouse retina expresses Zicl, Zic2, and Zic3 but expression of 
Zicl and Zic2 gradually disappears during retinal develop- 
ment, while Zic3 expression continues into adulthood 
(Watabe et al. 2011). Evidence that Zic factors are important 
for retinal development comes from overexpression of Zic 
genes in retinal progenitors, which decreases the number of 
rod photoreceptors (Watabe et al. 2011); thus, the increase in 
z/c7 and z/c2 expression in cavefish eyes may also contribute 
to the observed decrease in photoreceptors. 

Another category of genes that were downregulated in 
S. anophthalmus cavefish eyes involve catabolism, energy 
metabolism and ATP synthesis. This suggests that the degen- 
erated eye of Sinocyclocheilus cavefish has greatly reduced 
metabolic rate compared with surface fish eyes. Given that 
a functioning eye uses energy at a rate that even exceeds that 
of the brain (Wangsa-Wirawan and Linsenmeier 2003; Protas 
et al. 2007), selection may act to reduce this energy sink in 
species inhabiting environments in which retinal function is 
unnecessary for survival (Protas et al. 2007; Jeffery 2009). 

Down regulation of crx and Retinal Reduction in 
Cavefish 

Several studies have helped to improve our understanding of 
gene regulatory networks that influence photoreceptor de- 
velopment and homeostasis (Hennig et al. 2008; Hu et al. 
2010; Swaroop et al. 2010). RNA-seq data can test whether 
the regulation of genes in these networks was modified 
during cavefish evolution. The transcription factor crx is 
expressed in photoreceptors and other retinal neurons, and 
appears to be near the top of the hierarchy of retinal gene 
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regulation. Expression of crx enhances the expression of pho- 
toreceptor-specific genes and events downstream of otx2 
(Furukawa et al. 1997; Nishida et al. 2003; Shen and 
Raymond 2004; Swaroop et al. 2010) (fig. 7) and both crxa 
and its duplicate from the teleost genome duplication (TGD) 
crxb(otxS) are down regulated in cavefish eyes (fig. 6; supple- 
mentary fig. S3, Supplementary Material online). Knockdown 
of crxa in zebrafish results in eyes that are 10% smaller than 
control eyes and reduces expression of several phototrans- 
duction genes, including rhodopsin, cone opsins, irbp, otxS, 
gnatl, gnat2, and gucala (Shen and Raymond 2004). In Crx 
knockout mice, photoreceptors developed but failed to 
express many phototransduction genes, resulting in the 
deficiency of photoreceptor outer segments and eventually, 
retinal degeneration (Furukawa et al. 1999). Bipolar cells also 
express Crx, and their differentiation and maintenance is 
impaired in Crx knockout mice (Furukawa et al. 1999; 
Hennig et al. 2008). This evidence indicates that crx is near 
the top of the pathway for photoreceptor development and 
differentiation (Hennig et al. 2008; Swaroop et al. 2010). Like 
crx knockdown zebrafish and mice, Sinocyclocheilus cavefish 
also exhibit small eyes, short photoreceptor cells, and reduced 
numbers of retinal cells (figs. 1C and 2E-J). These features, 
combined with low crx expression and down regulation of 
genes downstream of crx (fig. 7; supplementary table S2, 
Supplementary Material online), are consistent with the 
mechanistic hypothesis that much of the histological eye 
phenotype of Sinocyclocheilus cavefish lies downstream of 
evolved decrease of crx activity that leads to shorter and 
fewer photoreceptor cells and reduced retinal cell gene ex- 
pression. Experiments to critically test the hypothesis that crx 
plays a direct role in retinal degeneration of cavefish requires 
further investigation. 

What causes the down regulation of crx expression in 
Sinocyclocheilus cavefish eyes? Although Otx2 is the only tran- 
scription factor yet known to directly interact with the pro- 
moter of crx and to regulate crx expression (Furukawa et al. 
1997; Nishida et al. 2003), otx2 expression was unchanged in 
the eyes of cavefish relative to surface fish. One hypothesis 
to explain these data is that cis evolution in the promotors 
of crx genes decreased responsiveness to Otx2 or other reg- 
ulatory proteins. Arguing against this hypothesis is that both 
crxa and crxb(otxS) were downregulated in cavefish eyes, 
and so mutations with similar effects would have had to 
occur in promotors of multiple paralogs, which is not parsi- 
monious. An alternative hypothesis is that cavefish have an 
evolutionarily mutated version of an as yet unidentified reg- 
ulator of crx genes, perhaps an inhibitor, that is responsible 
for the down regulation of crx paralogs in the cavefish 
S. anophthalmus. In ongoing research, we are looking for fac- 
tors that might regulate crxa expression during retinal devel- 
opment and maintenance. 

Comparing Sinocyclocheilus and Astyanax Cavefish 
Like Sinocyclocheilus, adult Astyanax cavefish have degener- 
ated eyes that are covered with skin and have sunken into 
their orbits (Jeffery et al. 2000; Yamamoto and Jeffery 2000). 



Also like Sinocyclocheilus, eye degeneration in A. mexicanus 
cavefish evolved several times independently in different cave 
populations (Tobler et al. 2010; Strecker et al. 2012). In 
Astyanax, retinal degeneration is initiated by apoptosis of 
the lens, which forms during embryogenesis and then degen- 
erates (Yamamoto et al. 2004). Our adult S. anophthalmus 
cavefish also have internal eyes (fig. 1C), which preserve basic 
but abnormal eye structures, including neural retina, iris, 
cornea, and lens (fig. 2B). In contrast to Astyanax, however, 
we did not find obvious abnormalities in either the morphol- 
ogy or transcriptomics of the lens in S. anophthalmus. Instead, 
the major morphological changes in Sinocyclocheilus cavefish 
eyes occurred in the neural retina. Although lens dysfunction 
is the causative factor in Astyanax cavefish eye degeneration 
(Yamamoto and Jeffery 2000), we find no detectable lens 
degeneration in Sinocyclocheilus cavefish. 

For transcriptomics, the Sinocyclocheilus transcriptome 
was annotated primarily based on the zebrafish genome 
while Astyanax transcriptomics used cross-hybridization 
to a zebrafish microarray (Strickler and Jeffery 2009). 
Phylogenically, Sinocyclocheilus (Cyprinidae) and zebrafish 
(Cyprinidae) are more closely related to each other than 
either species is to Astyanax (Characidae) (supplementary 
fig. S1, Supplementary Material online). Although some 
genes may have been missed in the cross-species microarray 
for Astyanax, for genes that did hybridize, comparisons 
between surface and cave dwelling Astyanax will be valid. 
Comparisons of gene expression profiles of cavefish/surface 
fish pairs for Astyanax (Strickler and Jeffery 2009) and 
Sinocyclocheilus revealed only four genes that were differen- 
tially expressed in both genera of cavefish, including cish, 
rhodopsin, gnatl, and gnat2. Cish (cytokine inducible SH2- 
containing protein), is an anti-inflammatory molecule that 
is upregulated in the cornea after interleukin stimulation 
(Ueta et al. 2011); in cavefish eyes, it is likely responding 
to inflammation associated with retinal dysfunction. The 
visual pigment rhodopsin and the transducin alpha subunits 
of rods and cones, Gnatl and Gnat2, respectively, are mem- 
bers of the visual phototransduction cascade (Shen and 
Raymond 2004). These comparisons suggest that retinal 
degeneration in A. mexicanus cavefish and S. anophthalmus 
share downstream effectors in retinal degeneration, although 
their reduced expression in S. anophthalmus might be 
caused by a lens-independent process, such as regulation by 
crx or other factors. 

Comparative analysis of morphology and gene expression 
patterns in two evolutionarily independent cavefish line- 
ages suggest that eye reduction in these two taxa appears 
to be controlled by different mechanisms, a lens-dependent 
process in Astyanax but a lens-independent process in 
Sinocyclocheilus. We do not know, however, in either case, 
the upstream evolutionary genetic changes or the selective 
regime, if any, that cause the retinal degeneration. 

Conclusions 

In this study, we contrasted histology and gene expression in 
the whole eye transcriptome of two closely related species: 
the surface species S. angustiporus and the cave species 
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S. anophthalmias. Results showed that cavefish developed 
distinct troglomorphic characters in the evolutionary adap- 
tation to perpetual darkness, including eye reduction, duck- 
bill-like snout and long barbels and pectoral fins. Cavefish 
eyes showed reduction in eye size, retinal cell density, photo- 
receptor cell height, and cell proliferation, and downregula- 
tion of genes encoding photoreceptor protein products 
and several retinal cell regulators, including crx, although 
expression of otx2, the only known direct crx upstream 
regulator, was unchanged in cavefish eyes. These findings 
suggest the hypothesis that the evolution of crx dysregulation 
by a novel trans-acting factor is a major contributor to the 
down regulation of retina-specific genes and reduced eye size 
in Sinocyclocheilus cavefish. 

Our analyses of Sinocyclocheilus cavefish retinas were lim- 
ited to samples from adult time points. As such, although we 
did not detect pyknotic nuclei or the presence of apoptotic 
markers in retinal tissue, we cannot exclude the possibility 
that cell death at an earlier time point contributes to the 
overall reduction in retinal cells observed. Our results, how- 
ever, do show substantial reductin in cell proliferation in adult 
retinas compared with surface fish, presenting the possibility 
that there are multiple mechanisms affecting the retinal 
changes we observe in Synocyclocheilus cavefish. Future inves- 
tigations will focus on determining which of the observed 
changes in gene regulation might influence reduced prolifer- 
ation or retinal degeneration, and at which developmental 
stage these phenomena initiate. Finally, the retinal reduction 
observed in Sinocyclocheilus cavefish appears to be indepen- 
dent of lens degeneration, in contrast to Astyanax cavefish. 
These results show that evolution can tinker with different 
developmental mechanisms to achieve similar adaptations to 
cave environments. 

Materials and Methods 

Animals 

Sinocyclocheilus fish were collected in Yunnan province, 
China. Cavefish (S. anophthalmus) were collected in Jiuxiang 
cave (in perpetual darkness), Yiliang County, and surface fish 
(S. angustiporus) were collected from Huangnihe River in 
Agang Town, Luoping County (supplementary fig. S5, 
Supplementary Material online). Both collecting sites are in 
the Nanpanjiang River drainage, the largest tributary of 
Xijiang River of the Pearl River basin. Although the two 
collecting sites are only approximately 30 km apart, they 
are isolated by a mountain range (supplementary fig. S5, 
Supplementary Material online). Cavefish were maintained 
in the laboratory in a dark environment and surface fish 
were exposed to a natural daylight cycle. All experimental 
procedures involving animals were conducted and approved 
by the Animal Care and Use Committee of Institute of 
Zoology, Chinese Academy of Sciences. 

Histological Analysis and Immunohistochemistry 
For H&E staining and immunohistochemistry, adult S. angu- 
stiporus and S. anophthalmus fish were euthanized, eyes were 
removed, cleaned of associated skin and adipose tissue, fixed 



in 4% paraformaldehyde in phosphate-buffered saline (PBS), 
and equilibrated in 30% sucrose prior to cryosectioning. To 
improve fixative infiltration to the retina, the cornea was 
perforated with a small pin. Slides were stained with H&E 
using standard procedures. For immunohistochemistry, slide 
edges were lined with a Gnome PAP Pen liquid blocker and 
washed three times for 5 min with PBS before blocking for 1 h 
in blocking solution (10% normal goat serum and 2% bovine 
serum albumen in 1 x PBS). We used the monoclonal anti- 
bodies Zpr1 and Zpr3 (Zebrafish International Resource 
Center, Eugene, Oregon, USA) (Vihtelic et al. 1999) to label 
unknown epitopes found in red/green double cones inner 
segments and the outer segments of rods and red/green 
double cones, respectively; antibodies against protein kinase 
C (Millipore, 05-983) to label ON bipolar cells in the INL; and 
antibodies against glutamine synthetase (Millipore, MAB302) 
to label Muller glial cells in the INL (Shkumatava et al. 2004). 
Antibodies were diluted 1:200 with blocking solution and 
incubated overnight at 4°C Slides were then washed three 
times for 15 min with phosphate buffered saline with Tween 
20 at 25°C We then applied anti-mouse Alexa Fluor-488 
(Invitrogen) secondary antibody solution (diluted 1:200 in 
blocking solution) to slides and incubated at room tempera- 
ture for 2-3 h in a dark humid chamber. Slides were coun- 
terstained with ToPro3 (Invitrogen, Cat no.: 642/661) to label 
DNA in cell nuclei, mounted in Vectashield (Vector, H-1000, 
stored at 4°C), and coverslips were sealed with clear nail 
polish. Fluorescent images were analyzed using a Bio-Rad 
Radiance 2100MP confocal microscope system. 

For analysis of proliferation and apoptosis, primary antisera 
were rabbit antibodies (Cell Signaling Technology, Inc.), as 
follows: for proliferation, anti-Phospho-Histone H3 (no. 
3377); for apoptosis, anti-Cleaved Caspase-3 (no. 9661), 
which used at the dilution of 1:1,000 and 1:400, respectively. 
Zpr3 labeled the photoreceptors. The secondary antibodies, 
Alexa Fluor-488-coupled anti-rabbit IgG and Alexa Fluor-594- 
coupled anti-mouse IgG (Invitrogen), were used at 1:200. 
Slides were mounted in Vectashield mounting medium 
with DAPI (Vector, H-1200). 

Quantification of nuclei and photoreceptor length mea- 
surements were performed with ImageJ 1.43u software. 
Nuclear densities and cell length were compared between 
surface and cave species using a t-test. All data met the 
assumption of normality. 

Constructing the Reference Transcriptome 
Total RNA was isolated using TRIzol reagent (Invitrogen) 
from whole eye and brain of adult surface and cave species, 
respectively. Two individuals were used for each species. 
Oligo(dT) selection was performed using MicroPoly(A) 
PuristTM (Ambion) according to the manufacturer's proto- 
col. Superscript III First-Strand Synthesis System (Invitrogen) 
was used to create first-strand cDNA. Second-strand cDNA 
was synthesized using random primers and Klenow exo-po- 
lymerase (Epicentre). cDNA was sheared to short fragments 
with Biorupter Sonicator, to which sequencing adapters were 
ligated. After size selection, suitable fragments were used as 
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templates for PCR amplification, which, although it might 
occasionally introduce some artifacts, is necessary to obtain 
sufficient quantities for analysis. Four libraries were prepared, 
brain and eye libraries for each of two species. Finally, 
two brain cDNA libraries of both species (200-400 bp 
insert) were sequenced using 80-nt paired-end reads from 
the lllumina sequencer GA-II, and two libraries of eye samples 
(200-400 bp insert) using 100-nt paired end reads on the 
HiSeq 2000 (lllumina Inc., San Diego, CA, USA). 

De novo assembly of the Sinocyclocheilus transcriptome 
was carried out using the short read assembly program Trinity 
(http://trinityrnaseq.sourceforge.net/) (Grabherr et al. 2011). 
After removing adapter sequences and low quality base calls, 
raw reads of surface fish and cavefish were assembled using 
default parameters with 300 nt minimum contig length. We 
made separate contig files of surface and cave species. 

To generate a common and nonredundant Sinocyclochei- 
lus transcriptome database for further RNA-seq analyses, we 
merged contig files from the two species using Cap3 (Huang 
and Madan 1999). If contigs were more than 95% identical 
and overlapped by more than 80 bp, they were merged. The 
longest contig of the group was retained in the transcriptome 
database, and redundant contigs were removed. Libraries of 
merged contigs were defined as nonredundant transcrip- 
tomes of Sinocyclocheilus for this article. 

Contigs were aligned by BLASTX to the NCBI Nr protein 
database using an E-value cut-off of 10 -6 (E value 
< 0.000001). Alignment results were used to retrieve gene 
names from genelaccession provided by NCBI. Annotation of 
the Sinocyclocheilus contigs was accomplished using the 
Blast2GO program by mapping the BLAST results to the 
GO database and finally selecting GO annotation using 
Blast2GO rules with default values (http://www.blast2go. 
com/) (Conesa et al. 2005). We also mapped the GO anno- 
tations to the Generic GO-Slim terms using Blast2GO. This 
Transcriptome Shotgun Assembly project for S. angustiporus 
has been deposited at DDBJ/EMBL/GenBank under the 
accession GAHO00000000. The version described in this arti- 
cle is the first version, GAHO01 000000. This Transcriptome 
Shotgun Assembly project for Sinocyclocheilus anophthalmus 
has been deposited at DDBJ/EMBL/GenBank under the 
accession GAHL00000000. The version described in this 
paper is the first version, GAHL01 000000. 

Differential Gene Expression Using RNA-seq Analysis 
After removing barcodes and low quality base calls, we 
mapped trimmed lllumina reads from surface fish and cave- 
fish eyes to our nonredundant Sinocyclocheilus transcriptome 
using Bowtie (Langmead 2010). Because of sequence variation 
between paralogs and orthologs, the maximum difference in a 
read was set to three nucleotides, which is the maximum 
number of differences supported by Bowtie; other parameters 
were default. 

Relative to some other cypriniforms, including zebrafish, 
the lineage leading to Sinocyclocheilus, to goldfish (Carassius 
auratus auratus), and to carps (including C carpio) experi- 
enced a whole-genome duplication (Ohno et al. 1968; Leggatt 



and Iwama 2003; Mable et al. 2011; Wang et al. 2012; Braasch 
and Postlethwait 2013) that we call the CGD. The lineage 
leading to teleosts, including zebrafish, which is the closest 
relative to Sinocyclocheilus with a sequenced genome, 
experienced an earlier genome duplication event, the TGD 
(Amores et al. 1998; Postlethwait et al. 2000; Taylor et al. 2003; 
Jaillon et al. 2004). In annotating the Sinocyclocheilus eye 
transcriptome, we collapsed the two copies from the CGD 
into one annotation unit corresponding to the zebrafish 
ortholog. 

To assess variation between Sinocylocheilus paralogs 
and their zebrafish orthologs, we randomly picked 16 genes 
that were abundantly expressed in Sinocyclocheilus eyes and 
analyzed the average identities of paralogs or orthologs 
between surface and cave species. Results showed that 
the average identity of paralogs duplicated in the CGD 
within Sinocyclocheillus species was 97.6 ±0.96%, and the 
average identity of orthologs between the two species was 
98.1 ± 0.90% (supplementary fig. S6, Supplementary Material 
online). Current high-throughput methods have difficulties 
distinguishing paralogs with such high identities from short- 
read RNA-seq experiments; thus, when counting gene expres- 
sion levels, we ignored the CGD paralogs in these two species 
so the relative expression value was the sum of the two CGD 
paralogs. In addition, identities between paralog-aa (or para- 
log~ab) and either of the two paralog-b genes {paralog-ba 
or paralog-bb) was about 81% in Sinocyclocheilus, similar 
to that of zebrafish paralog-a and paralog-b (supplementary 
fig. S6, Supplementary Material online), which were produced 
by the TGD (Amores et al. 1998). With these similarities, 
Sinocyclocheilus reads could be aligned to the correct paralogs 
at the level of the zebrafish genome. 

The average identity of all Sinocyclocheilus contigs with 
their best hits to zebrafish (D. rerio) was 90.5 ± 4.2% (supple- 
mentary fig. S6, Supplementary Material online). When we 
aligned transcriptome reads to the D. rerio transcriptome, 
the number of mapped reads was about half the number 
that mapped to the Sinocyclocheilus transcriptome. Because 
of evolution in the Sinocyclocheilus lineage after the CGD 
and to sequencing errors, we inferred that some valid 
Sinocyclocheilus sequence reads might not have been consol- 
idated into the de novo assembly. By BLASTing transcriptome 
contigs against the zebrafish UniGene database, we identified 
gaps in the Sinocyclocheilus transcriptome contigs. To 
account for this problem, we first mapped trimmed RNA- 
seq reads to the de novo Sinocyclocheilus transcriptome 
assembly using Bowtie; then, we rescued reads that failed to 
map to the Sinocyclocheilus transcriptome by aligning them 
to the zebrafish UniGene database. 

To assess differential expression, we merged the number 
of reads of surface fish or of cavefish that mapped to the 
same UniGene of either our Sinocyclocheilus transcriptome or 
the zebrafish UniGene database, and converted mapped read 
numbers (including reads mapped on both Sinocyclocheilus 
and D. rerio) to RPKM values (Mortazavi et al. 2008). 
To enhance statistical robustness, we excluded from analysis 
genes with fewer than five RPKM in either species, although 
these genes do appear in supplementary table S2, 
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Supplementary Material online. Differential expression analy- 
sis for RNA-seq data used a rigorous algorithm that computes 
a conditional probability for observing N1 reads for a gene in a 
treatment given that N2 reads were observed in controls and 
experimentals with no biological replicates (Audic and 
Claverie 1997). P values correspond to the significance of dif- 
ferential expression. Log 2 FC and P values were computed for 
each gene. If log 2 FC more than 1 and P value less than 0.05, 
genes were considered to be differentially expressed (Schartl 
et al. 2012). 

Real-Time PCR 

cDNA samples were created using the Superscript III First- 
Strand Synthesis System (Invitrogen) from the mRNAs of the 
eyes of different surface and cavefish adults than those used 
for the RNA-seq experiments. Real-time PCR was conducted 
using SYBR Green (TaKaRa) chemistry. Real-time PCR primers 
were designed based on locations with identical sequences of 
the relevant gene in both species derived from our de novo 
assembled Sinocyclocheilus transcriptome using PerlPrimer 
design software (http://perlprimer.sourceforge.net/) (supple- 
mentary table S4, Supplementary Material online). The beta- 
actinl gene was used as an internal control. All real-time PCR 
assays were run in a 96-well block with 20 jil reactions in each 
well. In each assay, reference samples were set up, and the 
expression levels of target genes were then quantified relative 
to this sample. A total of 16 genes were tested with Real-time 
PCR. Each reaction (samples and primers) was repeated at 
least three times. 

Phylogenetic Analysis 

In our transcriptomes, we identified nucleotide sequences of 
z/c 7 (z/c family member-1) and gprSS (C protein-coupled 
recepto-85 also known as sreb2), both of which are present 
in single copy in zebrafish and were previously used as nuclear 
phylogenetic markers in ray-finned fish (Li et al. 2007, 2008). 
For both Sinocyclocheilus species, two CGD paralogs of both 
genes were found and categorized as "a" and "b" paralogs 
based on shared nucleotide substitutions. For other ray- 
finned fish and for tetrapods, z/c7 and gprSS sequences 
were obtained from Genbank or Ensembl. See supplementary 
table S5, Supplementary Material online, for accession 
numbers. 

Sequence alignments were generated with MUSCLE (Edgar 
2004) for both genes individually and then concatenated. For 
Sinocyclocheilus, "a" paralogs (z/c 7 a + gprSSa) and "b" para- 
logs (z/c7b + gprSSb) were concatenated. A maximum likeli- 
hood phylogenetic tree based on 1,857 bp (z/c7: 903 bp; gpr8S: 
954 bp) was generated with PhyML (Guindon et al. 2010) 
using the GTR model and 100 bootstrap replications. 

Identification of Enriched Pathways and Diseases 
We used KOBAS 2.0 (Xie et al. 2011) to identify enriched 
pathways and diseases from among genes that were differen- 
tially expressed in cave and surface fishes. The sequences of 
differentially expressed genes were compared with the Homo 
sapiens database using the "annotate" feature in KOBAS 2.0 



to allow inferences from the abundant data on human path- 
ways and diseases. We then used "identify" to find signifi- 
cantly enriched pathways and disease-associated genes; 
"inputs" were the two "annotate" results of upregulated 
and downregulated gene sets and the "background" was 
the entire set of 9,649 unique genes expressed in 
Sinocyclocheilus eyes identified by RNA-seq. Data were statis- 
tically analyzed using the hypergeometric test and Benjamini- 
Hochberg FDR (false discovery rate) correction, and only 
pathways or diseases with corrected P value less than 0.05 
were considered to be enriched. 

Supplementary Material 

Supplementary tables S1-S5 and figures S1-S6 are available at 
Molecular Biology and Evolution online (http://www.mbe. 
oxfordjournals.org/). 
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